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ABSTRACT 

Modelling of the spectra of magnetic A and B main sequence stars is generally done assuming 
that all the lines are split by the magnetic field according to the Zeeman effect. However, a 
number of prominent spectral lines are produced by closely spaced doublets or triplets. Such 
r"| . lines should be treated using the theory of the partial Paschen-Back (PPB) effect. Depending 

Qli on the strength and orientation of magnetic field, the PPB effect can result in the Stokes / 

and V profiles of a spectral line that differ significantly from those predicted by the Zeeman 
effect theory. It is important to understand the size and types of errors that are introduced into 
magnetic spectrum synthesis by treating such lines with the usual Zeeman splitting theory 
rather than using the correct theoretical treatment of line splitting. To estimate the error intro- 
duced by the use of the Zeeman approximation, numerical simulations have been performed 
for spectral lines of the element silicon, for which a number of important lines are actually 
in the PPB regime, assuming an oblique rotator model, for various silicon abundances and 
V sin / values. A comparative analysis of the Stokes / and V profiles calculated assuming the 
PPB and Zeeman splitting has been carried out for a number of both strong and weak Si II 
and Si III lines. The analysis indicates that for high precision studies of some spectral lines 
the PPB approach should be used if the field strength at the magnetic poles is B v > 10 kG. In 
the case of the Si II line 5041 A, the difference between the two simulated profiles is caused 
by a significant contribution from a so called "ghost" line. The Stokes / and V profiles of this 
particular line simulated taking into account PPB splitting provide a significantly better fit to 
the observed profiles in the spectrum of the magnetic Ap star HD 318107 than the profiles 
calculated with Zeeman splitting. Employing the PPB approach, the Si II 5041 A line can be- 
come a useful tool for abundance mapping and reconstruction of magnetic field configuration, 
due to its sensitivity to the silicon abundance and to the magnetic field strength. 

Key words: atomic processes - magnetic fields - line: profiles - stars: chemically peculiar - 
stars: magnetic fields 



1 INTRODUCTION 

In the presence of a magnetic field, spectral lines become wider due 
to the splitting of the energy levels corresponding to the individ- 
ual transitions. The phenomenon was discovered by Pieter Zeeman 
( 18971, after whom it was named the Zeeman effect. Applying the 
experimentally determined splitting and polarisation characteristics 
of the Zeeman effect to observed solar spectra, George Ellery Hale 
( 1908 1 was able to demonstrate the presence of a magnetic field in 
sunspots. Using a similar technique for the analysis of stellar spec- 
tra, Horace Babcock dl947t discovered Zeeman splitting of absorp- 
tion lines in circularly polarized spectra of the A2p star 78 Vir and 
deduced the presence of a 1.5 kG magnetic field in its atmosphere. 
Nowadays, the Zeeman effect is widely used to estimate the mag- 



netic field strength and configuration in the stellar atmospheres of 
magnetic stars, some of which also show peculiar chemical abun- 
dances (Donati et al. 1997; Mathys & Hubrig 1997; Wade et al. 
2000a). 

When the Zeeman effect is observed in an atom in a mag- 
netic field that is small enough that the magnetic splitting of the 
energy levels responsible for the transition is small compared to the 
fine structure (spin-orbit) splitting of those levels, a wide variety of 
magnetic splitting patterns is observed. This behaviour is termed 
(for historical reasons) the "anomalous" Zeeman effect. However, 
when the field is strong enough that magnetic splitting becomes 
substantially larger than spin-orbit fine structure splitting, the ob- 
served magnetic pattern splitting changes gradually to a simple 
("normal") Zeeman triplet. The extreme case of magnetic splitting 
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in this large-field limit was first analysed by Friedrich Paschen and 
Ernst Back (1921 1, and was later named the Paschen-Back effect. 
Fifty years later Maltby ( 11971) discovered that the magnetic split- 
ting of Li I 6708A resonance doublet in sunspots shows a Paschen- 
Back pattern. 

In the analysis and modelling of spectral lines observed in 
magnetic upper main sequence stars, it is normally assumed that 
the splitting of the lines is correctly described by the anomalous 
Zeeman effect. However, in some cases the fine structure splitting 
of one or both levels involved in a transition is very small, and 
magnetic fields found in some stars are large enough to produce 
magnetic splitting comparable in size to this small fine structure 
splitting. In these cases, the usual anomalous Zeeman effect for 
isolated energy levels no longer describes the splitting correctly. 
Instead, the splitting of the line should be calculated taking into 
account both the fine structure splitting and the magnetic splitting 
simultaneously. This regime is known as the incomplete or partial 
Paschen-Back (PPB) effect. The partial Paschen-Back regime can 
occur when observed spectral line profiles are created by closely 
spaced doublets or triplets. A number of spectral lines which are 
important for analysis of magnetic upper main sequence stars have 
such structure, for example the triplet lines of He I, the 6707A res- 
onance lines of Li I, the 4481 A line of Mg II, and numerous lines 
such as XX 4128, 4200, 4621, 5466 and 7849 of Si II. In the pres- 
ence of a sufficiently strong magnetic field (typically some kG), 
the magnetic splitting for such lines can become comparable to or 
larger than the value of the fine structure splitting, and the resulting 
pattern for the split components can differ significantly from the 
one calculated according to the usual anomalous Zeeman effect. 

Thus, it is natural to suppose that examples of the partial 
Paschen-Back pattern could be found in the stellar spectra of stars 
with strong magnetic fields. Most such stars belong to the class of 
chemically peculiar (CP or Ap-Bp) stars, and for them it is impor- 
tant to correctly model magnetic line splitting in order to obtain 
realistic results from studies of magnetic field structure and from 
abundance analysis. The first study of how the Paschen-Back split- 
ting affects the Stokes profiles of spectral lines formed under the 
typical conditions of magnetic Ap and Bp stars was carried out by 
Mathys (1990) in his study of the line profiles of the Fe II multiplet 
74 doublet XX 6147-49 in stars with fields of several kG. A general 
discussion of the theoretical methods used to describe the transition 
from the Zeeman to the Paschen-Back effect in conditions appro- 
priate to such stars was presented by Landolfi et al. ( 2001 ). 

During the last decade, several articles have been published 
on the simulation of line profiles taking into account the partial 
Paschen-Back effect. Such work includes analysis of the Li 1 6708A 
resonance doublet QTochukhov 2008), of the He I 10830A multi- 
plet ( Socas-Navarro et al. 2004 ), and the components of the Fe II 
multiplet 74 (Stift et al. 2008). The influence of the Paschen-Back 
effect on the scattering polarization of molecular lines was consid- 
ered by Shapiro et al. (2006). 

Because of the importance of silicon in the spectrum of mag- 
netic CP stars, and the fact that many of the commonly analysed 
spectral lines of Si II are closely spaced fine structure doublets, we 
have carried out a study of the PPB magnetic splitting and of the 
resulting Stokes / and V profiles for a number of strong and weak 
Si II and Si III lines that are observed in the spectra of magnetic CP 
stars in the spectral band between 3 800 A and 8000A. Typically, 
such stars show enhanced silicon abundance in stellar atmospheres 
(Landstreet et al. 1989), and the location of their overabundance 
spots depends on the configuration of the magnetic field struc- 
ture (Hatzes 1991, 1997; Khochukhov 2002; Shavrina et al. 2010; 



Liiftinger et al. 2010). Therefore, these profiles are modified by the 
presence of a complex magnetic field structure and a non-uniform 
abundance distribution of silicon. Practical attempts to simultane- 
ously fit several silicon lines in the spectrum of the magnetic CP 
star HD 318107, which has a field of order 15 kG, show that, as- 
suming Zeeman splitting, some lines can be fit well in the frame- 
work of a chosen model for magnetic field structure, while the other 
lines, such as Si II 5041 A and 5056A, can not ( |Bailey et al. 201 1| >. 
This is an example of a modelling problem which suggests that the 
theoretical description of line splitting in a magnetic field could 
introduce unnecessary ambiguity into modelling. In an effort to un- 
derstand to what extent the difference between Zeeman splitting 
theory and that of the PPB theory is an important source of mod- 
elling errors, we have carried out a series of calculations to see how 
the Stokes / and V profiles differ when these lines are computed 
according to the partial Paschen-Back effect or the Zeeman effect. 

A description of the method for evaluating the difference be- 
tween simulated profiles assuming the partial Paschen-Back effect 
and the Zeeman effect is given in Sec. [2] and the results of its 
application to a number of Sill and Si III lines are presented in 
Sec. [3] Section|4]is devoted to analysis of the characteristics of the 
Sill 5041 A line magnetically widened in the regime of the partial 
Paschen-Back effect. The best fit results for the Stokes IV profiles 
of this line in HD 318107 spectra are also presented in that section 
for simulations that employ the partial Paschen-Back effect and the 
Zeeman effect. Discussion of the results obtained is given in Sec. [5] 

2 METHOD 

The theory behind the calculation of partial Paschen-Back splitting 
can be found in Chapter 3 of the recent monograph on polarization 
in spectral lines by Landi DeglTnnocenti and Landolfi (2004). The 
energy levels are calculated by diagonalization of the total Hamil- 
tonian, which consists of unperturbed and magnetic parts, under 
the assumption that the atomic system is described by the L-S cou- 
pling scheme. For the purposes of this paper, the unperturbed en- 
ergy levels (for zero magnetic field) of different ions are taken from 
the NIST atomic database (Ralchenko et al. 201 1 ) and the oscilla- 
tor strengths of the respective transitions are taken from the NIST 
and VALD (Kup ka et al. 2000[ > atomic databases. 

For weak magnetic fields the results of the Zeeman and 
Paschen-Back splitting computations are almost the same. With an 
increase of the field strength the magnetic splitting becomes com- 
parable to or even larger than the fine-structure separation between 
J-levels (where J stands for the total angular momentum of the 
atomic system). This condition gives rise to the partial Paschen- 
Back effect, which results in wavelengths and intensities of the 
split components which differ significantly from the predictions of 
the Zeeman effect for the same strength and direction of the ap- 
plied magnetic field. Furthermore, in the PPB regime transitions 
with |A7| > 1 are possible. In strong magnetic fields such forbid- 
den transitions produce spectral lines of relatively strong intensity. 
Such forbidden multiplet components can be observed in the spec- 
tra of some Ap stars (Mat hys 1990} and have been called "ghost" 
lines (Landi DeglTnnocenti & Landolfi 2004; Stift et al. 2008). 

In this work the numerical calculation of Paschen-Back split- 
ting is carried out using equation (3.64) from the aforementioned 
monograph. The code (by V.K.) has been designed for calculation 
of the Paschen-Back splitting and has been verified through com- 
parison of the results obtained with numerous examples (splitting 
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of the terms 2 P, 2 D, 2 F, 3 P, 3 d£] from Landi Degl'Innocenti and 
Landolfi J2004t . with the results of splitting of the Li I 6708A res- 
onance doublet ( Kochukhov 2008 ) and with previous computations 
of the components of the Fe II multiplet 74 ( Sti ft et al. 2008b . The 
splitting and relative intensities of lines from Table[TJin the regime 
of the partial Paschen-Back effect are consistent with the results ob- 
tained by using the code designed by Landi Degl'Innocenti (private 
communication) for calculation of Paschen-Back splitting. Thus, 
our method gives results similar to those previously published. 

2.1 Calculation of Paschen-Back splitting 

After the verification of the part of the code designed to calculate 
Paschen-Back splitting, it was incorporated into the ZEEMAN2 
code. This code was created by Landstreet ( 19881 for the simulation 
of polarimetric (Stokes IV QU ) line profiles, and was later modified 
by Khalack & Wade (2006), who added an automatic minimiza- 
tion of the model parameters using the downhill simplex method 
( IPressetal. 19921 ). 

To find the partial Paschen-Back splitting for even a small 
range of wavelengths, all the components of the multiplet under 
consideration must be taken into account because the intensity nor- 
malization of the split components is done over the whole multi- 
plet. The procedure for calculation of Paschen-Back splitting takes 
into account the unperturbed (or magnetically perturbed) energy 
levels and determines the respective air wavelength and oscillator 
strength of components, based on the term configurations and the 
total strength of all lines in the multiplet under consideration. With 
this code it is possible to simulate a line profile composed of several 
spectral lines (blends), some of which are split in the Paschen-Back 
regime, while the others are split in the Zeeman regime. One of the 
two routines (Paschen-Back or Zeeman) is applied to each of the 
contributions to the blends. 

To find the difference between the spectral profiles simulated 
assuming the PPB and Zeeman regimes, it is essential to require 
that both sets of computations are carried out with identical atomic 
data. Under the assumption of zero magnetic field there should be 
no difference between these profiles. The procedure for calculat- 
ing partial Paschen-Back splitting recalculates the wavelengths and 
intensities of all the line components for each field strength, includ- 
ing zero field. Therefore, for a correct comparison of the two cases 
it is necessary to adopt the relative line strengths within the mul- 
tiplet, as derived in the PPB case for zero magnetic field, for the 
calculation of Zeeman splitting of the same set of spectral lines. 
For the condition of zero magnetic field there is no actual split- 
ting in the PPB regime, but the relative intensity of each line in 
a multiplet is determined through the sum of its components as- 
suming the Russell-Saunders (or L-S) coupling scheme for the en- 
ergy levels (Landi Degl'Innocenti & Landolfi 2004). This same set 
of zero-field relative intensities must be adopted for the Zeeman 
splitting case. However, the total intensity of the multiplet is nor- 
malised to the sum of oscillator strengths of its members taken from 
the NIST or, sometimes, VALD atomic databases. The adopted os- 
cillator strength of each line is determined by multiplication of its 
relative strength in the multiplet by this sum (see the third column 
in Table [TJ(. 

The oscillator strengths derived from L-S coupling for spectral 
lines in the multiplets studied here generally show good agreement 

1 The splitting of the higher terms (up to 8 H) has also been calculated to 
verify the code and the similarity of splitting patterns. 



with the NIST data. However, the VALD data have been used for 
overall multiplet normalisation when they provided better agree- 
ment with the L-S coupling relative intensities of spectral lines in 
the multiplet studied. In Table[TJthe oscillator strengths from VALD 
are employed for simulation of Sill lines Xk 5041, 5056, 5957, 
5979, 6347 and 6371, and the remaining oscillator strengths are 
derived from the NIST database. Note that the perturbed part of 
the analyzed Hamiltonian includes only the contribution from the 
magnetic field and is nil when the field vanishes. All other possible 
contributions to the perturbed part of the Hamiltonian are not con- 
sidered here. Due to the simplified approach of the description of 
the analyzed Hamiltonian, some lines (Xk 3856, 3862, 4128, etc.) 
still show small differences between the values tabulated in NIST 
or VALD and the computed values of oscillator strength. 

To ensure that the "ghost" lines play a negligible role in the 
case of zero magnetic field, their oscillator strengths are set to have 
logg/ =-9.999 when the field is zero for the purposes of simula- 
tions. It should be pointed out that not all multiplets have "ghost" 
lines in the framework of the Paschen-Back splitting theory, be- 
cause some multiplets do not involve any pairs of energy levels 
with |A7| ^ 2. 

2.2 Comparison of the Paschen-Back and Zeeman profiles 

The spectral lines studied are listed in Table [TJ This table collects 
(separately for Si II and Si III) the multiplet number according to 
Adelman et al. ( 1985 1, the air wavelengths of the lines, the logg/ 
values obtained from the NIST and VALD databases (columns 4 
and 5) and the adopted oscillator strengths from L-S coupling the- 
ory, normalised to the best matching database values, in column 
3. The remaining columns list the basic atomic parameters of the 
lower and upper energy levels involved in each transition, and the 
Lande splitting factors used for Zeeman splitting computations. 

The simulation of Stokes / and V profiles is carried out for 
the lines in Table [TJ assuming a star with effective temperature 
r e ff=13000K, gravity logg=4.0, zero microturbulent velocity, and 
an oblique rotator model with a dipolar magnetic field structure. To 
explore the dependence of the final results on the field orientation, 
the calculations have been repeated for three cases where the axis of 
the magnetic dipole forms an angle oc= 0°, 45° and 90° with the line 
of sight, and the results are presented by the continuous, dashed and 
dotted lines, respectively (see Fig.[TJl. To investigate the influence 
of the axial stellar rotation and other non-magnetic phenomena that 
contribute to the widening of a line profile, the analysis is also per- 
formed for different values of Vs'mi (see Fig.(3](. Because most of 
the magnetic CP stars have strong Si II lines due to an enhanced sil- 
icon abundance (see, for example, Landstreet et al. 1989), the value 
log [Ns[/Nu] = —3.5 is adopted here for the simulations. The sili- 
con abundance is assumed to be homogeneously distributed in the 
stellar atmosphere. 

The observed Stokes / and V profiles are usually contaminated 
by observational noise. For practical reasons, it is useful to compare 
the difference between the simulated PPB and Zeeman profiles to 
the noise level O ~ at the analyzed spectral line. To decide 
whether this difference can be confidently detected above the noise 
level we use the chi square probability function and calculate the 
X 2 test statistics given by: 

where if and lf e represent the intensity of the Stokes / (or VQU) 
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Table 1. The air wavelengths and oscillator strengths of Si II and Si III lines analyzed, assuming no magnetic field. 
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profiles at a wavelength point i calculated with the assumption 
of the PPB and Zeeman splitting, respectively. The sum is eval- 
uated over the wavelength range within which the intensities of 
both profiles are significantly different from the continuum (or 
from zero in the case of Stokes VQU profiles), which means that 
\lf B \ > o and/or |/? e | > a. For a particular line this test quanti- 
fies how much its PPB profile deviates from its Zeeman profile 
distorted by the noise. If this distortion reaches the level of differ- 



ence between the PPB and Zeeman profiles the noise will mask it. 
The probability q that the simulated PPB profile is indistinguishable 
from the simulated Zeeman profiles plus noise with the dispersion 
O can be found from the chi-square probability function p(% 2 | v ) 
(Abramowitz & Stegun 1972) as: 

?(X 2 |v) = l- P (X 2 |v) = eQ,f ) , (2) 
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where v is the number of resolved elements in the analyzed profile 
and Q denotes the upper incomplete gamma function: 

For a small magnetic field the % 2 value is small and the proba- 
bility q that the noise masks the difference between the PPB and 
Zeeman profiles is high. Nevertheless, % 2 grows with increased 
field strength and for y} = 35.7 (assuming v = 20) the aforemen- 
tioned probability drops below 0.01. This confidence level seems to 
be quite robust to perform an evaluation of the difference between 
the PPB and Zeeman profiles for individual lines (see, for instance, 
Donati et al. ( |1992| >). If this criterion is satisfied (q ^ 0.01), the 
choice of magnetic splitting regime (PPB or Zeeman) during the 
simulation of the spectral line under consideration influences the 
results of spectral analysis. 

Since the ESPaDOnS spectropolarimeter is one of the most 
popular instruments designed to perform spectral observations of 
Stokes IV QU parameters, its spectral resolution R = 65000 and 
signal to noise ratio in polarimetric mod^| are adopted here for our 
analysis. The given spectral resolution provides approximately 20 
resolved elements in the Stokes / profile of the analyzed lines. This 
number increases with the magnetic field strength as the Stokes / 
profile becomes wider. 

Depending on the purpose of the spectropolarimetric obser- 
vations, the S/N ratio in the Stokes / spectra routinely obtained 
with the ESPaDOnS, NARVAL and MuSiCoS (the prototype of 
ESPaDOnS) is between 100 and 1000 (Folsom et al. 2010; Wade 
et al. 2000a). Using the ESPaDOnS spectropolarimeter some ob- 
servers even work with S/N~1300 (or more) in the Stokes / 
spectra in order to detect weak magnetic fields in massive stars 
(Grunhut et al. 2012) or in order to reach a higher precision in re- 
construction of a magnetic field configuration and/or abundance 
maps. To provide the same approach to the investigation of dif- 
ferent spectral lines, an average value S/N=250 (0=0.004) in the 
Stokes / spectrum is adopted here for all spectral regions. 



3 RESULTS OF ANALYSIS 

The list of Si II and Si III lines that are strongest in Stokes / spectra 
is given in Table[]]for a star with r eff =13000K and logg=4.0. The 
same list of lines is also used to analyze the Stokes V line profiles. 

A magnetically widened profile sometimes blends together all 
lines of a multiplet (for example Sill 5466A, 7849A, etc.). This 
demands their simultaneous analysis. In other cases the Stokes / 
and V profiles are analyzed separately for each line of the multiplet. 
In this study, the weak lines and the "ghost" lines are analyzed only 
if they contribute to a profile formed by a strong spectral line. 

3.1 General results for a slowly rotating star 

The difference between simulated profiles assuming the PPB and 
Zeeman splitting is first analyzed for a slowly rotating star assum- 
ing Vsin(= 1 km s and a silicon abundance of log [Ngi/Nn] = 
—3.5. The oblique rotator model with a dipolar magnetic field is 
applied to describe the configuration of a magnetic field. The pa- 
rameter a defines the angle between the dipole magnetic axis and 

For more details about this instrument, visit 
http://www.cfht.hawaii.edu/Instruments/Spectroscopy/Espadons/ 



the line of sight. The field strength at the poles of the magnetic 
dipole considered here changes gradually from zero up to 30 kG 
(for B p < 10 kG the step in the field strength is 500 G and for B p > 
10 kG it is 1 kG). Note that a polar field strength B p of 10 kG cor- 
responds typically to an observed mean field modulus (\B\) of 7 or 
8 kG, and to a mean longitudinal field strength (S-) of about 3 kG 
or less, so this polar field strength is typical of quite a number of 
known magnetic ApBp stars. 

Depending on the sensitivity of the probability logarithm (see 
Eq. |2j to the field strength at the magnetic poles, all the analyzed 
lines can be divided into two groups. The first group includes spec- 
tral lines for which the probability q calculated for the Stokes / pro- 
files decreases strongly with the magnetic field strength (see Fig.[TJl. 
Some of these lines are comparatively strong (Sill XX 3853, 3856, 
3862, 4128, 4130, 5041, 5056) with a relative intensity in Stokes / 
spectra of less than 0.6 (for V sin i= 1 km s ), while the others are 
weak (Sill XX 4200, 4621, 5466, 7849). Except for the Sill lines 
XX 3853, 3856, 3862, 4130 and 5056 all of them show an average 
difference between the Paschen-Back and Zeeman profiles that can 
be masked by the noise with the probability q < 0.01 (assuming 
S/N=250) in the Stokes / spectrum for a magnetic field strength 
S p =5 - 15 kG. For the Sill lines XX 3856, 3862 and 4130 there is 
no dependence of logg on magnetic field strength for Bp <20 kG 
(see the left panel of Fig.[TJl. 

The slope of the \ogq with respect to the strength of magnetic 
field is smaller for Stokes V profiles than for Stokes / profiles for 
these lines. Assuming the S/N=250 for the data in Stokes V spec- 
tra, the Si II lines XX 4128 and 5041 show an average difference be- 
tween the Paschen-Back and Zeeman profiles that exceeds the level 
of significance for Bp > 15 kG and Bp >7 kG, respectively (see the 
right panel of Fig. [TJ. In the case of Si II line 504lA the differ- 
ence appears mainly due to the presence of a comparatively strong 
"ghost" line in the PPB splitting (see Sec.[4jl. For the Stokes VQU 
spectra the S/N ratio averaged over the available line profiles is 
usually much less than the value obtained for the Stokes / spectrum 
(Donati et al. 1997). According to the definition of the Stokes pa- 
rameters (see for details chapter 1 in Landi DeglTnnocenti and Lan- 
dolfi ((2004)) the Stokes VQU spectra are obtained as a difference of 
two measurements taken with different combinations of orientation 
of a quarter-wave retarder and/or a polarizer (Donati et al. 1999). 
In this case the resulting signal is significantly smaller comparing 
to the Stokes / spectrum, while the average uncertainty remains the 
same. Therefore, if we adopt S/N=50 (CT=0.02) for the Stokes VQU 
spectra there will be no significant difference between the simu- 
lated Paschen-Back and Zeeman profiles (the probability q is close 
to 1). 

For weak magnetic fields there is no significant difference 
in the behaviour of the \ogq for different angles a. However, for 
Bp > 5 kG the \ogq decreases faster with the magnetic field 
strength for oc= 0° than for oc= 45° for both Stokes / and V profiles. 
The main reason is that the contribution to line profiles from the 
area of visible stellar disc that is close to the equator of magnetic 
dipole (where the field is weaker than near the poles), grows with 
the angle a producing higher % 2 (see Eq.QJ. For the case when oc= 
90°, there is no data for the Stokes V profiles (see the right panel of 
Fig. QJ, because the signal becomes too weak (smaller than a) for 
this orientation of the magnetic dipole. For the Stokes / profiles and 
a= 90° one can see that the decrease of the log q with the magnetic 
field strength is less steep compared the other two cases. 

The dependence of the equivalent widths (EW) on the mag- 
netic field strength for most of the lines listed in Table Q] is shown 
in Fig. [2] In general, for lines from the first group, the equiva- 



© 2010 RAS, MNRAS OOO.mmi 



6 V. Khalack & J. D. Landstreet 



Stokes I 



Stokes V 




10 20 30 10 20 30 



O 



-5 - 



O 



o 



O 




O 



O 



t — i — i — i — I — i — i — i — i — I — i — i — i — i — I — r 



Sill 3856A 



-- Sill 3862A 



I I I I I I I I I | I I I I | I 



> Sill 4128A -- Sill 4130A 



i i i i | i i i i | ' i i 'i i | i 



1 1 1 1 1 1 1 1 n 1 1 1 i+h 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



Sill 5466A 



I I I I | I I I I | I I h-|h- 



Silll 5739A 



J I I I I I I I I I I I I I I U I I L 



i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — r 



i i i i I i i i i I i i i i I i 



i i i i | i i i i | i i i i | i 




I I I I III I I | I I I l-|-4 





Sill 7849A 



I I I I I I I I I I I I I I I I 



-- Sill 6347A 



i . ■ . . i . ■ ■ . i 



10 20 30 10 20 30 



Bp, kG 



Bp, kG 



Bp, kG 



Bp, kG 



Figure 1. Logarithm of the probability that the noise with a = 0.04 (S/N=250) is masking the difference between the Stokes / and V profiles calculated 
assuming PPB and Zeeman splitting, assuming log [Nsi/Nji] = —3.5 and V sin; = 1 km s -1 . The continuous, dashed, and dotted lines correspond to the cases 
where the axis of magnetic dipole forms an angle a= 0°, 45°, and 90° with the line of sight, respectively. The horizontal dash-dotted line corresponds to the 
probability q=0.Ql. The bottom panels present the examples of Si II and Si III lines that show q ss 1 for both Stokes / and V profiles with little or no dependence 
on the field strength. 



lent widths of Stokes / profiles calculated employing PPB split- 
ting grow faster with the magnetic field (continuous line) than the 
equivalent widths of profiles calculated assuming Zeeman splitting 
(dotted line). The Si II lines XX 4128, 5041 and 5466 show a differ- 
ence between the equivalent widths of the PPB and Zeeman profiles 
AEW > 0.020 A for a magnetic field B p > 10 kG. This fact indi- 
cates that when the magnetic field at the magnetic poles reaches 
a strength of 10 kG, the difference between the PPB and Zeeman 
profiles in unpolarized spectra for the aforementioned lines is quite 
significant (see also Fig. |7J- 



The spectral lines for which the logq shows very weak or no 
dependence on the magnetic field strength for the Stokes / and 
V profiles (see the bottom panels of Fig. QJ compose the second 
group. One reason for the occurrence of such spectral lines is that 
for transitions between terms with spin quantum number S=0, the 
pattern of the PPB splitting is the same as the one obtained in 
the Zeeman splitting {Landi Degl'Innocenti & La ndolfi 2004}. The 
Sim lines 5739A and 4716A (see Tabl.Q}, for which % z = for 
all values of the magnetic field strength (see Eq. Q]and the corre- 
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Figure 2. The equivalent widths of the analyzed lines as a function of magnetic field, assuming a silicon abundance of log [Nsi/Nu] = —3.5. The continuous 
and dotted lines correspond to the profiles calculated using PPB and Zeeman splitting respectively. For Sim lines 4716A and 5739A note that the vertical 
scale is ten times smaller. 



sponding panels for Si III line 5739A in Fig.QJ are good examples 
of this rule. 

The other lines (Sill U 5957, 5979, 6347, 6371 and Si III 
XX 4552, 4567, 4574) provide logq close to zero with almost no de- 
pendence relative to the magnetic field strength for the Stokes / and 
V profiles (see panels for Si II line 6347A in Fig. Q}. These lines 
belong to widely spaced multiplets, and so even for field strengths 
B„ ~ 30 kG they are still well in the Zeeman regime. We see clearly 
that spectral lines belonging to the second group can be analyzed 
without any loss of information assuming Zeeman splitting even 
for stars with very strong magnetic fields. 

The Stokes Q and U profiles are quite sensitive to the con- 
figuration of the magnetic field and to the abundance inhomo- 
geneities ( Wade et al. 2000a). Nevertheless, with the assumption of 
a S/N=250 for Stokes Q and U spectra and a homogeneous distri- 
bution of silicon abundance, the Stokes Q and U profiles appear to 
be rather weak. The analyzed silicon lines do not show any signif- 
icant differences between the Stokes Q and U profiles computed 
with the Zeeman effect or the PPB effect for B v < 30 kG and they 



can be simulated assuming the Zeeman splitting without any loss 
of information. With the assumption of a S/N<50 the Stokes Q and 
U profiles are barely detectable, hence in further analysis, most at- 
tention is focussed only on the Stokes / and V profiles computed 
according to the Zeeman and PPB effects. 

3.2 Dependence on V sin i 

The majority of magnetic CP stars show velocities of axial rotation 
with Vsim greater than 1 km s , and mostly less than 50 km s 
(Kudryavtsev et al. 2007; Wade et al. 2000b). The dependence of 
the magnetic field strength at which \ogq = —2 on the projected 
rotation velocity is studied here for some lines in order to evaluate 
its influence on the % 2 (see Eq. [TJ. The strength of the magnetic 
field determined in this way provides the threshold above which 
the line profile simulated in the PPB regime becomes significantly 
different from the one simulated assuming Zeeman splitting. 

The values of the polar magnetic field corresponding to 
logq = —2 are shown in Fig.[3]fbr several lines assuming S/N=250. 
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Figure 4. Difference between the Stokes / profiles calculated with the PPB and Zeeman splitting for the Si II lines 5041 A (left panel) and 5466A (right panel) 
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of the magnetic dipole. The dashed horizontal lines present the level of significance for p=0.99 (q=0.01, see Eq.[2). 
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Figure 3. Variation of the threshold of the magnetic field (see details in 
Subsec. |3~2l with V sin i for the Stokes / spectra. The open circles, triangles, 
and squares correspond to the cases where the axis of magnetic dipole forms 
an angle a= 0°, 45°, and 90° with the line of sight, respectively. 



For the Sill lines XX 4128 and 5041, the derived threshold grows 
slowly with V sin i for the Stokes / spectra and corresponds to field 
strength around 10-18 kG and around 5-10 kG, respectively. For 
Si II lines XX 4621, 5466 and 7849 the derived threshold increases 
as a function of Vsim, but remains less than 18 kG for Vsini < 
15 km s _1 . Meanwhile, for Sill 4200A the derived threshold in- 
creases rapidly as a function of V sin i. Similar results are derived 
for two other values of the angle a (see Fig.[3}- 

The results obtained suggest that it may be important to take 



into account the PPB effect during the simulation of the Si II lines 
of the first group, especially the line 5041 A, for magnetic stars with 
V sini < 15 km s~' 



andfi p > 10-18 kG. 



3.3 Dependence on silicon abundance 

Some magnetic CP stars show patches of an enhanced silicon 
abundance in their stellar atmospheres (Hatzes 1997; Khochukhov 
2002; Shavrina et al. 2010; Luftinger et al. 2010; Bailey et al. 
2011). An increase of silicon abundance leads to stronger line pro- 
files when keeping the other parameters constant. Two examples 
of how the dependence of logg on the magnetic field strength 
changes with increasing silicon abundance (homogeneously spread 
over stellar surface) are shown in Fig. [4] for the Sill lines 5041 A 
(left panel) and 5466A (right panel), belonging to the first group 
(see Subsec. mt . For clarity only the curves that correspond to the 
model with an angle a = 45° between the line of sight and the 
axis of the magnetic dipole are shown, while the values of adopted 
silicon abundance log [N^/Nu] are placed on the panels near each 
curve. Similar behaviour of the logg is found for other values of 
the angle a. 

In the case of Sill 5041 A line, an increase of silicon abun- 
dance leads to a steeper dependence of the logarithm of proba- 
bility logq on the magnetic field strength and correspondingly to 
a smaller magnetic field at which logq = —2 and the difference 
between the PPB and Zeeman profiles becomes significant. This 
behaviour can be explained as the effect of intensification of the 
"ghost" line (see Sec. [4] and Fig. [7} with increasing silicon abun- 
dance, which in turn results in a bigger difference between the PPB 
and Zeeman profiles (larger % 2 in Eq. Q}. A similar behaviour can 
also be seen in the case of the weak line Si II 5466A. Neverthe- 
less, for a silicon abundance log [Nsi/Na] = —4.0 the dependence 
of the log q on the magnetic field strength changes its slope for B p > 
20 kG, while for log [Nsi/Nn] = —4.5 the logarithm of probability 
log g-O for Bp < 30 kG. 

The first example illustrates the unique behaviour of the Si II 
504lA line, while the second example shows a typical example of 
how a dependence of logq on the magnetic field strength changes 
with increased silicon abundance. The results obtained here indi- 
cate that variation of silicon abundance by about 1 dex does not 



© 2010 RAS, MNRAS 000,[TJtT2] 



Partial Paschen-Back splitting of Sill and SMI lines 9 



eliminate the difference between the PPB and Zeeman profiles; in 
general, this difference grows with increased silicon abundance. 



4 THE SI II LINE 504lA 

For most of the Si II and Si III spectral line profiles studied, here 
a contribution from their respective "ghost" components remains 
relatively small when the magnetic field strength is B p < 30 kG. 
The same result was found by Stift et al. ( 2008 1 during the analysis 
of PPB splitting of the Fe II multiplet 74. 

However, this is not true for the Sill line 5041 A. From the 
left panel of Fig. [4] we can see that the PPB effect becomes im- 
portant for simulations of the Stokes / profile even for solar sili- 
con abundance when the field strength at the pole exceeds 14 kG. 
Fig. H] shows the wavelength splitting and relative intensities of the 
subcomponents of the Sill line 5041.0A and the "ghost" line at 
5040.7A. It appears that at B p ~ 12 kG those components start 
overlapping, while the Jl and a lec |-components of the "ghost" line 
increase in intensity rapidly with the strength of the magnetic field. 

For a smaller magnetic field (fi p < 12 kG) the main differ- 
ence between the Paschen-Back and Zeeman profiles comes from 
the depression caused by the "ghost" line. As an example, Fig. [7] 
shows the comparison of Paschen-Back and Zeeman profiles in the 
Stokes / (left panel) and V (right panel) spectra simulated assum- 
ing enhanced silicon abundance log [Ns[/Nu] = -3.5 with B p = 10 
kG and rx=45° for most of the analyzed lines. In the case of Sill 
504lA line the "ghost" component is located around V= - 20 km 
s and its contribution to the left wing of the line is significant. 
It increases with the silicon abundance (see left panel of Fig. |4j 
and can be incorrectly identified as a blend if one uses the Zee- 
man effect approximation to simulate the Stokes / and V profiles 
of this particular line. A spectral analysis of the depression caused 
by the "ghost" component is a good tool for estimating the sili- 
con abundance and determining the magnetic field configuration of 
stars with spectra clearly showing the presence of silicon lines. 

4.1 The 5041A line in the spectrum of HD 318107 

To test the theoretical results obtained for the Sill line 5041 A the 
Stokes / and V profiles of this particular line and of the line Si II 
6347A were analyzed for the magnetic Ap star HD318107 using 
both Zeeman and PPB splitting. The 6347A line is chosen as a 
comparison, because its Zeeman pattern of split components does 
not differ much from the pattern of PPB splitting for any magnetic 
field strength B p < 30 kG, and correspondingly the parameter d is 
close to zero (see Sec. 13. II ). The polarimetric spectra were obtained 
at the rotational phase 0.991 using the echelle spectrapolarimeter 
ESPaDOnS at CFHT with a resolving power of R ~ 65000 and 
signal-to-noise ratio S/N ~ 200 (Bail ey et al. 201 1\ . 

HD318107 is a very peculiar magnetic star with a rotational 
period of P = 9.7088 ± 0.0007 days, deduced by Bailey et al. 
( 1201 11 1 from analysis of the periodic variability of its magnetic 
field components. For the rotational phase analysed here its mean 
longitudinal magnetic field (S ; ) is around 5 kG and the mean 
field modulus (|S|) is about 14 kG. This star has an effective 
temperature T eS = 11800±500 K and gravity logg= 4.22±0.13 
dLandstreet et al. 2 007 1. The stellar atmosphere model was calcu- 
lated for temperature r e ff= 11800 K and gravity logg= 4.2 using 
the PHOENIX code dHauschildt et al. 19971 assuming LTE and so- 
lar metallicity. The modified ZEEMAN2 code (see Sec. I2.lt was 
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Figure 5. (Colour online) Splitting and relative intensities of the subcom- 
ponents of the Sill lines 5041.0A and 5040.7A ("ghost" line) in the PPB 
regime. The magnetic field strength of 10 kG is represented by the vertical 
line and the Oblue, 1, and le j components are plotted respectively in blue, 
green, and red. 



employed to analyze the spectropolarimetric data assuming a ho- 
mogeneous distribution of silicon abundance and a magnetic field 
geometry characterized by a sum of colinear dipole, quadrupole, 
and octupole at an angle (3 to the rotational axis. To describe the 
geometry of a magnetic field we have adopted from Bailey et al. 
( f20TTl the following parameters : (=22°, p=65°, B d = 25.6 kG, B q = 
-12.8 kG, B = 0.9 kG. 

Results of the best fit assuming Zeeman magnetic splitting for 
both lines and PPB splitting for the Sill line 5041 A are shown 
in the Table [2] where the chi-squared values are given for com- 
bined analysis of the Stokes / and V profiles. For each case listed 
in this table the fitting of silicon abundance, V sin(/) and radial ve- 
locity has been performed individually using the downhill simplex 
method (Pr ess et al. 1992t . In the case of the Si II 6347A line, ap- 
plication of the Zeeman splitting approximation results in a good 
fit (see Fig. [6j, while in the case of Sill 5041 A line, PPB split- 
ting provides a much better fit than can be obtained with a Zeeman 
profile. 

The weak blending lines Fel 5040.85A, 5040.90A and Fell 
5040.76A also can make a contribution to the observed profile in 
the area of the "ghost" component of Si II 504lA line. To analyze 
their influence on the best fit results, these lines were included in the 
computations, and split in the Zeeman regime during the simulation 
routine, while the Si II line 504lA was split in the PPB regime. 
The respective best fit results are given at the second line of Table|2] 
(see dash-dotted line at the Fig. [6jl and are much better in the PPB 
regime than the results obtained for the same set of lines split in the 
Zeeman regime (dotted line). 

One can see that the data obtained independently for both lines 
agree between themselves and are close to the results derived by 
Bailey et al. d201 lb from a complex analysis of different spectral 
lines for this star. Nevertheless, the depression in the left wing of 
the Sill 5041 A line can not be fit sufficiently well taking into ac- 
count its PPB splitting and the contribution of the iron blends. The 
remaining differences between the PPB profiles and the observed 
spectra may be partially explained in terms of more complicated 
magnetic field structure and/or horizontal and vertical stratification 
of the silicon abundance ('Bailey ~et al. 201 1) . 
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Figure 6. The best fit model for the Stokes / (left) and Stokes V (right) profiles of Sill 5041 A and 6347A lines assuming Zeeman (dotted line) and PPB 
(dashed line) splitting. Continuous lines represent the observed data, while the dash-dotted lines provides the best fit results for Si II 504lA in the PPB regime 
in combination with the Fel and Fell blends in the Zeeman regime for log [Npe/Nji] = —3.0 . In the case of Sill 5041 A the dotted line shows the best fit 
results for the silicon and iron lines split in the Zeeman regime. The results for the Zeeman splitting of this line without contribution from the iron lines are not 
shown here. The 6347A spectra are shifted vertically for better visibility. 



Table 2. Approximation of Sill 5041 A and 6347A lines observed in the 
spectrum of HD 318107 (phase 0.991) employing the PPB and Zeeman 
effects for their magnetic splitting and adopting from Bailey et al. t201H 
the following parameters for the magnetic field geometry: £=22°, p=65°, 
B d = 25.6 kG, B q = -12.8 kG, B = 0.9 kG. 



X, 


Splitting 


logpVsi/Ak] Vsin(i), 


Vr, 


X 2 /v 


A 




km s -1 


km s _1 





5041 


PPB 


-3.57i0.15 


7.9±1.0 


-8.4=1=1.0 


16.5 


5041 


PPB+Fe 


-3.80±0.15 


7.3=1=1.5 


-8.1±1.0 


9.8 


5041 


Zeeman 


-3.74=1=0. 15 


7.5±0.8 


-8.0=1=1.0 


70.7 


5041 


Zeeman+Fe 


-3.94±0.15 


7.5±0.8 


-7.7=1=1.0 


40.4 


6347 


Zeeman 


-3.67=b0.15 


7.4±1.0 


-8.9=1=1.0 


7.8 



5 SUMMARY 

This paper presents the results of simulation of Stokes / and V pro- 
files of strong and weak Si II and Si III lines (see Tab. QJ in the 
presence of stellar magnetic field. The purpose of this study is to 
find out the conditions under which the application of the Zeeman 
effect approximation is unreliable. 

Presently, the Zeeman effect is used in most of the codes de- 
signed to calculate the Stokes IV QU profiles of spectral lines. For 
magnetic CP stars with a comparatively weak magnetic fields its 
application is generally valid for most of the spectral lines and re- 
sults in a correct reconstruction of the abundance maps and the 
magnetic field configuration. Even for the stars with comparatively 
strong magnetic fields (20 kG < B p < 30 kG), whose spectra have 
a S/N ratio around 50, the Stokes / profiles of almost all the Si II 
lines (except 4 128 A, 5041 A and 5466 A) may be treated with the 
Zeeman approximation. Nevertheless, if the available polarimetric 
spectra of the stars with strong magnetic field have higher S/N ratio 
(for example S/N=250 and more) the use of PPB splitting during 
the analysis of spectral lines is necessary to obtain precise results 
in the framework of an assumed model for the abundance map and 
the magnetic field structure. In particular, Stokes / profiles of the 
Si II lines which belong to the first group, as defined in Subsec. [JTl 
when calculated with the PPB splitting, differ significantly from 



those calculated with the Zeeman effect. This difference appears 
due to the different relative intensities and positions of split a and 
7t-components in the PPB and Zeeman regimes, and due to the so 
called "ghost" lines (\AJ\ 2) as in the case of Sill 5041 A line 
(see Sec. [4]!. 

In our simulations (except Subsec. 14. It we adopt a signal to 
noise ratio S/N= 250 in the Stokes / and V spectra. A portion of 
our analysis has been performed for an enhanced silicon abundance 
log [Nsi/Nn] = —3.5, Vsinz'= 1 km s . Increasing the rotational 
velocity causes a little decrease of the difference between Paschen- 
Back and Zeeman profiles for some Si II lines and diminishes it sig- 
nificantly for others. The analysis of the Sill lines Xk 4128, 5041, 
5466 and 7849 shows that the application of the PPB splitting for 
the simulation of Stokes / profiles remains important in stars with 
V sini < 15 km s when the field strength is B p > 10-18 kG. 

The weak silicon lines can usually be used for spectral analy- 
sis only if the spectra are obtained with a high S/N ratio or if silicon 
is significantly overabundant in the stellar atmosphere. For higher 
S/N ratio (S/N> 250), significant differences between the Paschen- 
Back and Zeeman profiles can also be found for weaker magnetic 
fields. The higher silicon abundance causes a steeper decrease of 
log q (the logarithm of probability that the noise masks the differ- 
ence between the PPB and Zeeman profiles, see Eq.[2]for details) 
with the magnetic field strength (see Subsec. 13.3b . If the patches 
of enhanced silicon abundance (approximately by 1 dex) provide a 
significant contribution to the observed line profiles, the PPB effect 
should be applied to correctly simulate the magnetic broadening 
of the Stokes IV profiles of the silicon lines belonging to the first 
group, even for a star with a moderate magnetic field (B p > 10 kG). 
This statement is especially important for the rotational phases of a 
star with a magnetic dipole axis pointing towards the observer. 

For the case of the the Sill 5041 A line most of the difference 
between the Zeeman and Paschen-Back profiles is caused by a sig- 
nificant contribution from the "ghost" line, which is well visible 
in the spectra with S/N=250 for B p > 6 kG and an enhanced sili- 
con abundance (see the right panel of Fig. [TJ. The contribution of 
this "ghost" line can grow essentially with the strength of magnetic 
field. If one relies on the Zeeman effect during the simulation of this 
spectral line, the observed depression in its left wing can only be 
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Figure 7. The simulated Stokes / (left panel) and V (right panel) profiles of some silicon lines, applying the PPB (continuous line) and Zeeman (dotted line) 
splitting for an abundance log [N$i/Nu] = —3.5 and a magnetic dipole with the field strength of 10 kG at the magnetic poles. The spectra are shifted vertically 
for better visibility by 0.3 and 0.1 in the left and the right panels respectively. In the left wing of Si II line 504lA a contribution from the "ghost" component 
is clearly visible. 



partially explained by a presence of iron blends (see Subsec. 14.11 ). 
In contrast, the use of the PPB effect approach allows for the in- 
corporation of the "ghost" line and results in the much better fit of 
the observed in HD 318107 Stokes IV profiles. A combination of 
the iron blend split in the Zeeman regime and the Si II 5041 A line 
split in the PPB regime provides an even better fit (see Fig. [6] and 
Table [2}. Therefore, application of the PPB splitting of this partic- 
ular line provides us with a good tool for estimation of the silicon 
abundance and for the reconstruction of the magnetic field config- 
uration. 

The Paschen-Back effect approximation has not been com- 
monly used for simulations of the Stokes IVQU profiles mainly 
because it requires more sophisticated calculations and a high pre- 
cision knowledge of the energy levels for the transitions studied. 
Its use is becoming more common in recent astrophysical research, 



especially for the investigation of stellar objects with strong mag- 
netic fields. Application of the Zeeman effect to calculate the mag- 
netic broadening of all spectral lines in these objects can lead to 
imprecisions in some cases. A combination of the PPB effect anal- 
ysis for the first group lines with the less demanding Zeeman effect 
analysis for the second group lines can provide more precise re- 
sults. For stars with somewhat weaker fields, the important result 
of this work is to quantify the level of error introduced by using the 
Zeeman approximation in situations where the PPB effect should 
in principle be used, making it possible for investigators to decide 
which method is required in specific applications. 

The use of this technique for the analysis of stellar spectra will 
be explored by the authors in the near future as well as the similar 
analysis of spectral lines produced by the other chemical species in 
which the PPB effect may play a significant role. Meanwhile, for 
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the silicon lines of the second group an application of the Zeeman 
effect provides the correct results of magnetic splitting for compar- 
atively strong magnetic field. 
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